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Navier-Stokes equations are known as hydrodynamic equations which take account of effects of 
dissipations. There are, however, problems in the relativistic Navier-Stokes equations, i.e. the 
equations violate causality. Israel-Stewart equations, which evade the problems of Navier-Stokes 
equations by introducing new parameters, such as the relaxation times, have recently been used 
in describing the space-time evolution of the quark-gluon plasma produced in high energy heavy 
ion collisions. The viscosities and the relaxation times are related to each other by imposing 
entropy constraints on the system. According to Boltzmann-Einstein principle, the probability 
distribution of the fluctuation in the energy-momentum tensor is related to the entropy of the 
system. Applying this principle to the entropy in Israel-Stewart theory, one can obtain the ratios 
of the viscosities to the relaxation times. We evaluate the ratios of the viscosities to the relaxation 
times in SU(3) lattice gauge theory. 
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1. Introduction 

The new form of matter created at RHIC (Relativistic Heavy Ion Collider) has provided physi- 
cists unexpected findings. One of the most impressive results is that the new matter, which is 
referred to as quark-gluon plasma (QGP), behaves like perfect fluid near critical temperature T c . 
The space-time evolution of QGP in high energy heavy ion collisions seems to be understood 
quantitatively by relativistic ideal hydrodynamics. However, relativistic dissipative hydrodynamics 
is needed for more quantitative description of the real QGP evolution since QGP has small but 
nonzero viscosities. 

The simplest relativistic dissipative hydrodynamics is the first order theory which is a rel- 
ativistic extension of Navier-Stokes theory [[T], g]. This theory, however, accompanies problems 
associated with causality violation [Q]. Instead of the first order theory, second order dissipative 
hydrodynamics, especially Israel-Stewart (IS) theory [Q], has been recently actively studied. IS 
theory avoids the problems by introducing second order terms of dissipation, which physically rep- 
resent relaxation effects to the solution of the first order theory. As an inevitable consequence of 
introducing the second order terms, however, IS theory includes many phenomenological parame- 
ters in addition to the three transport coefficients in the first order theory. These parameters cannot 
be determined within hydrodynamics. Other approaches based on microscopic theories are needed 
to determine the parameters. 

"Perfect fluid"-like behavior of QGP indicates that quarks and gluons strongly interact with 
each other in QGP near T c The properties of QGP near T c are not within the reach of per- 
turbation theory, since perturbative expansion is poorly converging at large couplings. At present, 
lattice gauge theory is the only systematic approach which can calculate physical quantities in such 
a non-perturbative region. There have been done several attempts to measure transport coefficients 
with numerical simulations on the lattice [0, [7L [8]|. These analyses evaluate transport coefficients 
using Kubo formulae, which relate transport coefficients to the correlation functions of the energy- 
momentum tensors in Minkowsiki space-time. On the other hand, on the lattice one can calculate 
correlation functions in Euclidean space-time. An analytic continuation from imaginary-time cor- 
relation functions to real-time ones therefore is required in this strategy. In previous studies [||, [7|, M, 
nontrivial ansatze have been adopted for the form of real-time spectral function in order to perform 
this analytic continuation. Validity of such ansatze, however, should be carefully examined. 

In the present study, we attempt to constrain phenomenological parameters in IS theory on the 
lattice using a method proposed in Refs. [Kj. In this method, ratios between the viscosities and 
the relaxation times in IS equations are related to static fluctuations of stress tensor in equilibrium 
through Boltzmann-Einstein principle. Since one does not use temporal correlation functions in 
this method, one can avoid the difficulty in the analytic continuation. The main objective of this 
work is to evaluate these ratios in SU(3) gauge theory on the lattice and to reduce the number of 
phenomenological parameters in IS theory. 

2. Israel- Stewart theory 

In this section, we overview relativistic dissipative hydrodynamics, especially Israel-Stewart 
(IS) theory [Q, pT| ]. Basic equations of hydrodynamics are local conservation laws of the energy- 
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momentum and the net charge, 

5^7^ = 0, (2.1) 
d M W = 0. (2.2) 

Using arbitrary 4- velocity of fluid normalized as u^ 1 u ii = \ and the projection onto 3-dimensional 
space A^ v = g^ v — u^u v , the energy-momentum tensor T^ v and the charge density are decom- 
posed as 

= euW - (p + n)A" v + W"u v + W V + K^, (2.3) 
N^ = nu^+V^, (2.4) 

where £, and « are the energy density, pressure, and charge density, respectively. Eckart used the 
particle flow as 4- velocity , 

vf = , (2.5) 

which is called Eckart (particle) frame, and vanishes in this frame. Landau-Lifshitz employed 
the energy flow for , 

«" = v =■ (2.6) 

This is Landau-Lifshitz (energy) frame, and W^ v vanishes in this frame. While the values of 
and depend on the choice of the frame, 

=W H -Z+Ly* (2.7) 
n 

i.e. heat flow in the particle frame, is frame independent in first order. 



The basic idea of Israel [ 12 ] for a phenomenological derivation of second order hydrodynamics 
is to incorporate the effects of dissipation into the entropy current . Under the hydrodynamic 
assumption, i.e. that nonequilibrium states are characterized by hydrodynamic variables, T^ v and 
, the entropy current should be given by = s^{T pa ,N P ). Assuming further that one can 
expand for nonequilibrium states by the power series of IT, %^ v , and q^ , the most general form 
of entropy current at second order in dissipative terms, IT, n^ y , and q^ , reads, 

^ = ^ + C + G M > ( 2 - 8 ) 

n T 

where 

Q>—£ (fttf - Am' +**»»■*) - + 0.9) 

represents the second order contribution to entropy, and s eq and T are the entropy density in equilib- 
rium and the local temperature of the system, respectively. Here, /3, and a, are phenomenological 
coefficients and are not determined by the hydrodynamic assumption. 

Requiring the second law of thermodynamics, > 0, one can constrain the macroscopic 
equations that hydrodynamic variables follow. If we neglect the second order term in s^ 1 , the 
second law and linearity lead to the first order hydrodynamic equations including three transport 
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coefficients, shear and bulk viscosities 77 and £, and heat conductivity A. Similarly, when one 
recovers the second order terms, Q?, in Eq. ( |2.8[ ), the second law leads to constraints including 
second order terms. IS equations are defined so as to satisfy these constraints [Q, |T^ ]. 

IS equations include terms DTI and Dn with D = u^d^ being the derivative along uP- . These 
terms give rise to relaxation effects to the solution of first order equations, and cure the causality 
problem of the first order theory. Proportional coefficients of these terms are given by 



(2.10) 



respectively, which physically represent the time scale of the relaxation in each channel and are 



called relaxation times. Equation (2.1C) means that the ratios between viscosities and the relaxation 
times are related to proportional coefficients in the entropy current, jSq and p2, as 



ft - Tn 

PO — - J 



(2.11) 



3. Boltzmann-Einstein principle 



In what follows, we try to evaluate the ratios in Eq. ( [2.1 1| ) on the lattice by connecting fto and 
p2 to fluctuations of stress tensor using Boltzmann-Einstein (BE) principle. 

The entropy and the number of microscopic states in equilibrium are related to each other by 
Boltznmann relation 

S{a)=\nW{a), (3.1) 

where a = {01,02,03, ...} represents the set of state variables of the system, and S(a) and W(a) 
are the entropy and the number of microscopic states, respectively, in the system specified by state 
variables a. To put it another way, the number of microscopic states is given by BE principle 



10], 



W(a) = exp[S(a)]. 



(3.2) 



Since the probability that one of the state variables a, takes a certain value, A, in equilibrium is 
given by 

W(oi=A) 



P{ai = A) 



(3.3) 



ZW(d) ' 

BE principle means that the probability distribution of state variables in equilibrium is governed by 
the entropy, provided that the entropy S(a) is defined. 



>From Eqs. ( |2.8| ) and (|2.9|), the entropy in unit volume s in IS theory at rest frame uP is given 



by 



UnS r 



?eq - — (ft n 2 - p iq ^ + v^ v ) 



(3.4) 



Notice that with orthogonality relations u^q^ = u^n^ v = the terms including a, in Eq. ( pL9| ) do 
not appear in Eq. (^4|). As a metter of course, the entropy in Eq. ( 3^4 ) is the same quantity as that 



appearing in the BE principle. Substituting Eq. ( p.4| ) into Eq. (p.2|), the number of microscopic 
states specified by II, q^, and %^ v in a volume V is given by 



V 



W(Tl,q^7tn - exp[-— (poll 2 - frq^ + M^ v )], 



(3.5) 
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where n, q^, and are spatial average of each quantity in V, 



n = - / ^xn(x), (3.6) 



v 



V 

and so forth. Equation (EO) means that fluctuations of n, qf 1 , and 7T^ V are related to j3, in entropy 



current Eq. (2.8). In particular, distribution of each quantity is of Gaussian form; for example, that 



of jty v is given by 



and the fluctuation reads 



P(*">exp[-^ : J3 2 ji„ v f v ], (3.7) 

<(*") 2 > = ^, 0-8) 

where we have used (7T^ V ) = 0. 

Since j3,- are related to fluctuations of stress tensor through Eq. (3.8), one can evaluate the r] to 
Tjt ratio on the lattice by measuring fluctuations ((7T^ V ) 2 ) in a quite straightforward way. Although 
this analysis cannot determine the values of transport coefficients themselves, the number of free 
parameters in IS theory can be reduced. Since this method does not refer to correlation functions, 
one can avoid the procedure of the analytic continuation, which is one of the most nontrivial task 
in the measurement of viscosities using Kubo formulae. 

4. Measurement of fluctuations on the lattice 

The basic idea of lattice gauge theory is to discretize the four dimensional Euclid space in a 
gauge invariant way. The partition function on the lattice in path-integral representation is given by 

dUdifrdYerSoW-SFtoYP) = J dU DetFp] e^^. (4.1) 

Here Sq and Sp are the Euclidean gauge and fermion actions on the lattice, respectively. U, y/, and 
Y correspond to the gauge, anti-fermion, and fermion fields in the continuum limit, respectively. 
DetF[£/] represents the fermion determinant. In a general procedure of lattice calculations, the 
expectation value of a physical quantity O in equilibrium is computed as an average over gauge 
configurations: 

(0) = -( dU DetF[U] e- Sa{u) 0. (4.2) 

In this work we analyze SU(3) pure gauge theory where the fermion part in Eqs. (4T) and (4.2) is 
neglected; 

Z = JdU e- s ^ u \ (°} = ^J dU z~ Sg[U) 0. (4.3) 



In this theory, conserved current corresponding to N 11 in Eq. ( PL2j ) does not exist and the hydro- 
dynamic equation is solely given by Eq. ( |2.1[ ). Difference between energy and particle flows, q^, 
cannot be defined in this theory. 

On the lattice, the energy flow does not exist in the rest frame of the lattice. For n and 7fy v 
we thus take n = ij^Li Tn and n l -> = T,j with 1 < i < j < 3. Spatial average in Eq. (3.6) should 
be taken for a volume V on a given time slice. To eliminate possible artificial effects arising from 
periodicity along spatial direction, V should be considerably smaller than the lattice volume in 
numerical simulations. 
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j3=6/s 2 


a [fin] 


Size (Nl x N?) 


T/T c 


N config 


case 1 


6.499 


0.049 


32 3 x6 


2.5 


10000 


case 2 


6.499 


0.049 


32 3 x 10 


1.5 


10000 


case 3 


6.872 


0.031 


48 3 x 10 


2.5 


10000 



Table 1: Simulation parameters. N a , N z , and A^ con fig represent the number of lattice sites in spatial and 
temporal directions, and gauge configurations, respectively, a and N T are related to the temperature T as 
aNx = T~ l . Cases 1 and 2 have the same lattice spacing with the different temperature. Cases 1 and 3 have 
the same temperature with the different lattice spacing. 
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Figure 1: Distribution of fluctuation of 7Z\2 obtained from numerical simulations for easel. The horizontal 
axis is normalized by the temperature. The dotted line denotes the Gaussian fit. 



5. Numerical Results 

In this work, we perform the lattice simulation for SU(3) pure gauge theory with the standard 
Wilson gauge action. Configurations are generated by heatbath and overrelaxation algorithms. In 
Table [j] we show lattice parameters used in this work. The simulation is performed with three 
lattice sizes x N t with a periodic boundary condition and different lattice spacing a. For each 
parameter, 10000 configurations have been prepared. The energy momentum tensor is given by 
7^ v = 2tr [Fj?FP v + (1 /4)g^ v FP a F pa ] . For the definition of the field strength on the lattice, we 
have chosen the clover operator. 

For the spatial average in Eq. (3.6), we take the average of T^ v on a cube of size (N a /2) 3 to 



remove the effect of periodicity. For each configuration one can take eight such subconfigurations 
on cubes. We regard two of the cubes aligned in a diagonal direction as statistically independent 
areas. Furthermore, we adopt two sets of subconfigurations on the two time slices at N T = and 
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Figure 2: Distribution of 7Ti2 2 obtained from numerical simulations for case 1. The vertical axis is logarith- 
mic, fa can be read from the gradient of fit line. 



N t = 1 /2aT in order to improve statistics. In this way, four subconfigurations are taken from each 
configuration for analysis. 

In Fig. [lj we show the distribution of the spatial average of the shear viscous pressure ftn = T'yi 
for case 1 . The same result is shown in Fig. [| with different choices of vertical and horizontal axes. 
These figures show that the distribution of %\i is of Gaussian form, as anticipated. Extracting 
fluctuations from this distribution and substituting the result to Eq. (j3~8|), we obtain 



fa ^ 0.0036. 



(5.1) 



This value of fa is, however, quite small in the sense that the relaxation time x n = Ifar] is so 
small for T] < 1 that IS theory turns back to the first order theory. In fact, the speed of transverse 



plane wave, vj, with the value of fa in Eq. (5.1) gives [11] 



1 



vt_ = / 

c ~\2fa{e + p) 



(5.2) 



where we have used the values of e and p obtained in the previous work on bulk thermodynamics 



of SU(3) gauge theory [13]. Eq. ( pJ\ ) means that causality is violated with fa obtained in this 
analysis. Similar results are obtained on cases 2 and 3. 



6. Discussions 



In this work, we tried to evaluate the ratio between the relaxation time x n and the shear viscos- 
ity 17 in SU(3) gauge theory by lattice simulations. Using BE principle and the form of entropy in IS 
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theory, the ratio is related to statistical fluctuation of stress tensor in equilibrium through Eq. (3.8). 
In the numerical simulations, however, we have observed large fluctuations of stress tensor, (#i 2 )> 
which leads to unacceptably small /?2 violating causality. 



It is, however, suspicious to simply apply fluctuations obtained on the lattice to Eq. ( |3.8[ ), since 
(tt^ 2 ) is actually an ultraviolet divergent quantity and it diverges even in the vacuum. The large 
fluctuations would have originated from this vacuum fluctuations, which should be subtracted to 
obtain the physical value of fluctuations. This analysis is in progress. The present numerical results, 
however, indicate that we need more statistics to separate physical fluctuations from the vacuum 
one with a good accuracy. 
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